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Abstract 

Constraints on the required flatness of the scalar potential V{(j)) for a cousin-model to 
extended inflation are studied. It is shown that, unlike earlier results, Induced-gravity 
Inflation can lead to successful inflation with a very simple lagrangian and A ~ 10~^, rather 
than 10"^^ as previously reported. A second order phase transition further enables this 
model to escape the 'big bubble' problem of extended inflation, while retaining the latter's 
motivations based on the low-energy effective lagrangians of supergravity, superstring, and 
Kaluza-Klein theories. 



PACS numbers: 98.80C, 04.50 



1 Introduction 

Since Guth's first paper on inflation a decade ago there has been an explosion of effort dedicated 
to developing a natural theory of inflationary cosmology. These efforts have established inflation 
as a ubiquitous program of research among cosmologists and particle physicists, but have not yet 
produced any particular model which satisfies all of the practitioners. As Kolb and Turner recently 
remarked, inflation is by now "a paradigm in search of a model." Q Much of the trouble concerns 
the so-called "flne-tuning" required to force the various models into agreement with measurements 
of the small anisotropy of the cosmic microwave background radiation. 

Adams, Freese, and Guth |Q developed a quantitative means of measuring the amount of "flne- 
tuning" necessary to make inflationary models agree with observations. They deflne a dimensionless 
parameter A as the ratio of the change in the scalar fleld's potential to the change in the scalar 
held: 

A - (1) 

In this way, A functions as a measure of the flatness of a given potential. (Ref. |^ denotes this ratio 
by A, rather than A; the upper-case letter is used here to avoid confusion with the closely-related 
quartic self-coupling parameter.) In they evaluate this ratio for generic inflation scenarios (with- 
out a curvature-coupled scalar field) and for the original version of extended inflation [Q]. They 
show that extended inflation requires a fine-tuning eight orders of magnitude more stringent than 
the general inflation schemes: whereas they flnd A < 10~^ — 10^^ for the inflaton potential of new 
inflation, they calculate that a potential for the Brans-Dicke-like scalar fleld in extended inflation 
would require A < 10~^^. Although the authors of Q are quick to point out that their method of 
calculating A is highly model-dependent, others have taken this result to indicate that adding a 
potential V{(j)) for the curvature-coupled scalar in any type of extended inflation model necessarily 



entails this sort of extreme fine-tuning (see, e.g., [^). 

Tliis paper examines the constraints on V{(j)) for a cousin-model to the original extended in- 
flation. Building on Zee's early ideas Q about uniting spontaneous symmetry breaking with the 
Brans-Dicke reformulation of general relativity, Accetta, Zoller, and Turner ||7| developed a model 
of "Induced-gravity Inflation." In this model, a single scalar boson does all the work of inflation: it 
couples to the scalar curvature R and drives inflation with V{(t)) (unlike ordinary extended inflation, 
which requires one field to couple to R while a separate and unrelated field drives the expansion). 
This kind of one-boson model can agree with observations with about the same degree of "fine- 
tuning" as the generic models of inflation examined in Q: as will be shown below, "Induced-gravity 
Inflation" requires A < 0(10"^) - a far cry from the 10^^^ of ref. [^! Furthermore, by employing a 
second order phase transition to exit the inflationary epoch, rather than the first order transition 
of the 'big bubble' or 'w problem' which plagued original extended inflation |^ [l^ may be 
avoided. 

In section 2, we calculate A for Induced-gravity Inflation, and consider why earlier attempts to 
determine the required flatness of the potential have led to much more constrained results. Section 
3 examines the accuracy of the slow-rollover solutions upon which the calculation of A is based. 
And in section 4, we briefly consider benefits and difficulties of placing the inflationary epoch at 
such a high energy scale. 
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2 Calculating A 

We begin with the lagrangian density]^ 

= fim - \g^''d^(t^d,cl) - V{cl)) + Lm, (2) 

where the Lm term here only includes contributions from 'ordinary' matter - there is no separate 
Higgs sector as in Q. If we choose /((/;) = cjP' /{^uj), we find the coupled field equations: 

#+3/f0 + ^ = ^^4,-^[(p-3p) + iV{4,)-4,V'm- (3) 

In eq. (3), a{t) is the cosmic scale factor of the Robertson- Walker metric, and is related to the 
Hubble parameter hy H = a/a. From eq. (2), one can see that /(</>) — > (IGvrGejj)"^, which leads to 
4u;0~^ = SnGeff- Note that we have parametrized the lagrangian slightly differently from the 1985 
"Induced-gravity" paper: our Brans-Dicke parameter lv is inversely proportional to their coupling 
strength e {to = (4e)^^). That paper also makes the minor approximation that (2a;)/(3 -|- 2uj) 1, 
whereas we have kept this term explicit. During the inflationary epoch, the k term becomes 
negligible; similarly, since p and p now only include contributions from ordinary matter, they too 
may be neglected. If we now impose the appropriate slow-rollover approximations 0, |8|: 

< H, 

« V{^), (4) 



^The sign conventions follow those of Misner, Thorne, and Wheeler \^M, which are based upon the 1962 edition 
of Landau and Lifshitz 0. Thus, gtt < 0, the full Riemann tensor is = d^T^^ - a„r^^ + T^^F"^^ - T^^F"^^, 

and the Ricci tensor is i?^^ = ^^ak- The original Brans-Dicke papers followed these sign conventions. Note that 
these definitions for the Riemann and Ricci tensors are opposite in sign from those in Weinberg's text (and thus 
opposite to some of the recent literature on extended inflation). 
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the field equations enter a more tractable form: 



H 



2 ^ ^^yi<t>) 



3 



h2 



3 + 2u} (j) 

In order to calculate the required flatness (A) of the potential, we must now assume a particular 
form for V{(j)). To achieve a second order phase transition with a minimum of fine-tuning, we may 
choose the simplest form for the potential at tree-level: 

V{ct^) = \{<t^\t)-v^)\ (6) 

where A is the quartic self-coupling, and v is the vacuum expectation value for (p. In addition to 
being a qualitatively simple form for the desired potential, this is also the optimal form found in 
for the case of generic new inflation. It also matches the potential of [0, with the minor difference 
that it is parametrized with A/4 rather than A/8. 

Using this expression for V{(j)), combined with the field equations (5), we may solve for (j){t) 
and a{t): 



(7) 



The factor 7 = {3 + 2uj)/2 is slightly different from the exponent found in Q (because of their earlier 
approximation that (2w)/(3 -|- 2uj) — > 1). The quantities (po and as are values at the beginning 
of the inflationary epoch. Note that at early times, when (p{t) ~ 0o, these equations yield the 
familiar power-law solution, with a{t) oc As t increases, the rate of expansion slows due to the 
exp (</>o — 0^) term. With these analytic expressions for (j){t) and a{t), we may now calculate A. 
Following H, we begin with the two basic constraints on the scalar potential: it must provide 



for sufficient inflation to solve the flatness, horizon, and monopole density problems, and it must 
allow for the proper amplitude of density perturbations to act as seeds for the evolution of large 
scale structure. The first of these requirements takes the form: 

Hn aN ~ Hb a_B ' ^ ^ 

where the subscript N refers to present values ('Now'), and B refers to values at the beginning of 
the inflationary epoch. Using the standard cosmological model's assumption of adiabatic expansion 
following the end of inflation, we may write: 

O'end Tj<4 

-— - ^> (9) 

where Ognd is the value of the scale factor at the end of inflation, and Trh is the reheat temperature 
following thermalization of the foregoing false vacuum energy density. (Ref. assumes Trh ~ Mp, 
where Mp is the false vacuum energy density of the second boson which drives inflation.) The 
Hubble parameter may then be parametrized as: 

rr2 _ 87T Vf 

B — ~5 2' 
6 



where Vf is the initial value of the false vacuum energy density; nip is the effective value of the Planck 
mass at the beginning of inflation (related to the initial value of the field 0: (po = {io/2ir)^^'^mp); 
and Mp is the present value of the Planck mass, Mp M 1.22 X 10^^ GeV. The quantity (3 is the 
ratio of the energy density in matter to the energy density in radiation today, which Q takes to 
be around 81. In the original extended inflation framework, V = Vj = Mp was constant, due to 
the energy density of the metastable state of the second boson before it completed its first order 
phase transition. In the present model, V((/>) changes in time as (j){t) changes; yet even in this new 



context, Vf is still a constant, and is simply equal to V{(j)o); its use here is thus independent of any 
slow-roll approximation. 

Using eqs. (8-10), the condition for sufficient inflation takes the form: 



> 



1/2 



/ 



as PTmTrh rup 



(11) 



Combining with eq. (7), we get: 



as 



■'end 



exp 



^end 



V, 



1/2 



> 



/ 



(12) 



where we have followed Q in utilizing the explicit analytic solution for a{t) based on the slow-roll 
approximation right up to Qend] as we shall see below, this approximation is a good one: the slow- 
rollover solutions employed above don't break down until 4){t) ~ 0.98u. 

Eq. (12) may be used to place a bound on the change in the scalar field, — {4'end — ^o")- 
Writing (po and (pend in the exponent in terms of their associated masses, we find: 



^end 



> 



1/7 



exp 



-LO 



{ml -Ml,) 



(13) 



where we have followed Q in writing Me for the value of the Planck mass at the end of inflation. 

In the simplest case, (j) would not evolve any more after the end of inflation, so Me would equal 

Mp; yet for the time being the more general value Me will be used. Since 0o = (aj/27r)^/^?n,p and 

V ~ (pend = (i^/27r)-^/^M^, eq. (13) may be rewritten: 

/ r ..1/2 _ _ 1 1/7 



> \l — TTlp 



\ 



PTnTrh rup 



exp 



Me 



(14) 



Taking AF = Vf and defining /i = {ijip/Me), the ratio A thus becomes: 



A < 



27r\2 Vf 



I3TnTrh rUp 



1/7 



exp 



(15) 
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Incorporating the constraint that Vf must exceed the kinetic energy of the (j) field's de Sitter space 
fluctuations without exceeding the effective Planck scale at the beginning of inflation Q , we may 
bound Vf by: 

V,<(^)'< (16) 
which may then be used to eliminate the coefficient of Vf/rrip-. 



A < 





2 


/ - 




) 








V - 



I5TnTrh rup 



1/7 



exp 



1 



(17) 



Thus far we have only relied upon the requirement of sufficient inflation. The second require- 
ment, based on the amplitude of density perturbations, may now be used to place a bound on the 
ratio Tfj^jj /vj^'^ . This bound can then be used to remove nearly all mass-scale dependence from A. 
The post-COBE parametrization of the density perturbation constraint may be written as [pT|]: 



^ ^ J \hor 



(18) 



where the quantity H'^/(j) is to be evaluated at the time of last horizon-crossing of the perturbations 
(during the inflationary epoch), and 6h — 1-7 x 10~^. (Ref. Q used 5 x 10"^, rather than birdH-) 
As noted in [12|, the quantities pertaining to the time of last horizon-crossing in eq. (18) should 
be evaluated in the Einstein frame; when compared with the "naive" calculations in the Jordan 
frame, a correction factor F(u}) must be introduced. Yet this factor F{uj) decreases monotonically 
with Lo, and converges to unity rather quickly: F{lo = 25) = 1.052, and F{uj = 500) = 1.002. For 
the particular values of uj with which we shall be concerned below, F{uj) would thus be negligible, 
and so we may continue to calculate H'^/cj) in the Jordan frame. (This is also the approach adopted 
in 1^.) Since we want to relate eq. (18) to a bound on Vf, we should first rewrite cl){t^) in terms of 
Vf, where t* is the time of horizon-crossing. From eq. (7) we have: 



37' 



to 



2n 



(19) 



and from eq. (6) we have: 



Vf = yiM = ^Ml: {l-fi'f. (20) 
Using these expressions, we may then rewrite (the time-independent) (p in terms of the constant 



/2 



(21) 



The absolute value for the (1 — ^^) term comes from taking the positive square root of M|. bmce 
nip < Me, the pole at = 1 is exchided. 

Now we need to calculate the Hubble parameter at the time t*. Following we may write: 



1 



(22) 



which, after employing eqs. (9, 10, 21), leads to: 



HQ 



^end 



Rewriting this as a bound on (T|^/l^y^), we get: 



rp2 



< 



yl/2 j - 4^2 ^2 Y ^2 



(1 - /u^) < ^ttSh- 



(23) 



(24) 



This may be substituted into eq. (17) for A: 



A < 



4uj 



AI3Tn Trh 7' 



Oiend 



1/7 



exp 



(25) 



We have taken the equality in eq. (24) as a worst case for A; if (t^j^/Vj^^^ were much less than the 
right-hand side of eq. (24), the bound on A would increase (and the resultant need for "fine-tuning" 
would therefore decrease). The remaining factor of Trh/itt-p may now be removed by combining 
eqs. (16, 24). Following these substitutions, the expression for A becomes: 

2^l/2 



'37^'^' 



10it5h 



O'end 



1/7 



exp 



5 0--^) 



(26) 



The time of last horizon-crossing can be calculated from eq. (22), and from one could then find 
the ratio aend/o-{t^). Yet the dependence of A on aend/o-{t^) is very weak (being suppressed by the 
exponent 7~^), so some simplifying approximations may be made. In most models of inflation, t^, 
is around 60 e-folds before the end of inflation (although some recent models have ~ 50 e-folds 
before the end of inflation, e.g. [^); this means that the ratio aend/cL{t^) is simply e^^ (or, perhaps, 
e^'^). For the calculation of A, we will assume aend/ci{t^) ~ e^^ in this model. As Figure 1 shows, A 
increases monotonically with increasing for a given value of a lowest bound on A thus comes 
from taking the limit /i — > {i.e., nip <C Me). When this is done, A becomes a function of lo alone. 
Figure 2 shows a plot of A versus ui in the limit /x ^ 0: the maximum value of A (corresponding 
to the least amount of "fine-tuning" required) is 5.4 x 10~^, for the value oJcr = 240. We can check 
the dependence of A on a^nd/O'ii*) by defining a parameter a as dend/o-it*) = e". Figure 3 shows 
a plot of A versus a in the limit fj, ^ for a particular value of w (w = 500). The dependence 
on a is indeed weak: A evaluated at (a = 50,u; = 500) gives 4.5 x 10~^, whereas A evaluated at 
(q = 60,u; = 500) gives 3.8 x lO^^. Similarly, for q = 50 rather than 60, A{ujcr) = 7.3 x lO^^, 
instead of 5.4 x 10~^. 

In the original model of extended inflation, oj was constrained to be less than 25 in order to avoid 
observable inhomogeneities coming from the large range in bubble sizes, even though present tests 
of Brans-Dicke gravitation versus general relativity limit to to the range lo > 500 (hence the 'big 
bubble' or 'u problem' of old extended inflation). Yet in the present model, the second order phase 
transition lifts this constraint on uj; uj can now be as large as necessary to meet the experimental 
limits. As one can see in Figure 2, A{uj) falls off slowly from its maximum with increasing lo. It 
is interesting to note that Amax occurs within a factor of 2 of the value lo = 500. Because of its 
agreement with present day observations, and its proximity to ujcr, ^ = 500 appears to be a good 
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candidate for the Brans-Dicke parameter. 

Eq. (1) may be used to relate A to the quartic self-couphng constant A. For the form of V{(j)) 
considered here, we find: 

A = ^ = ^M^ (27) 
or, if we keep terms only up to 0{fj.), 

A^l + 0(^ (28) 

Thus, A is of the same order of magnitude as A, which, for the model under study, means 
A ~ 0(10"^). This value of A is much larger than the results in 1^] for the original model of 
extended inflation, indicating far less of a need for "fine-tuning" . 

We should pause here to consider why this relatively large value for A has not been noted be- 
fore. The most important reason is because the calculation of A is highly model-dependent. Both 
papers of |^], for example, assumed that a separate Higgs sector would drive inflation; this meant 
that their vl^'^ = Mp was constrained to lie at the GUT scale, with such ratios as {Mp/v) ~ 10 . 
Furthermore, by insisting upon a first order phase transition in the Higgs sector, u) was constrained 
to cj < 25. It is interesting to note that an attempt in 1989 to unite the original Induced-gravity 



Infiation model with extended infiation |16] similarly relied upon a separate Higgs sector to drive 
inflation until it underwent a first order phase transition. 

The authors of the 1985 paper introducing Induced-gravity Inflation Q studied constraints on 
the quartic self-coupling A, based also on the twin requirements of sufficient infiation and a proper 
amplitude for density perturbations. Yet their result indicated that A < 10~^^ for oj ~ 500. Several 
factors help to explain this low result. First, their pre-COBE parametrization of the amplitude 
of density perturbations leads to an increase of an order of magnitude for A when compared with 
present, post-COBE values. Most important, however, is their approach to bounding A: they 
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solved for A in terms of the ratio where (/>(t*) is the value of the field at the time of last 

horizon-crossing. Their result for A, which in their analysis is proportional to sinh^^[ln{v / ipit^,))], 
is thus very sensitive to the value of (j){t^). Because this value cannot be solved exactly (even in the 
slow-roll approximation), two extreme limiting regimes were studied: (j){t^) ~ v versus (/>(t*) <C v. 
Yet small differences in the approximation of (j){tif) lead to order-of- magnitude differences in their 
estimation of A: a difference of 0.01 in the assumed value of (/>(t^,) leads to a difference in A of 
two orders of magnitude. The method of calculating bounds on A employed in this paper avoids 
expanding in terms of the unknown ratio {v / (j){t^:)); the only mass ratios involved here are of order 
{(po/v)'^, and their inclusion raises the bound on A. The information regarding (pit^) is now con- 
tained in the ratio (aend/o-i't*)), and we saw above that changing this ratio from e^^ to e^^ leads to 
a change in A by a factor of only 1.18 (see Figure 3). This appears to be the major reason for 
the large split in values of A between this paper and the 1985 analysis. 



3 Accuracy of Slow-Rollover Approximate Solutions 

We may now check the accuracy of our slow-rollover approximate solutions by following Steinhardt's 



and Turner's "prescription" for successful slow-rollover [15|. The analysis is easiest by rewriting 
eq. (2) in terms of a Brans-Dicke field where ^ = /(</>) = (jP'/{^uj). The "prescription" of |15| 



concerns finding conditions for when the $-term may be neglected. In the present model, when $ 
is negligible, the $ equation becomes: 
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where the prime now indicates differentiation with respect to Using this expression for we 



may calculate and then write the ratio $/(3i?$), which becomes: 



3H^ 9iJ2 (3 + 2io) 



^ [V - ^V"] - ^ 



1 f9H\ n 



(30) 



Prom eq. (30), it is consistent to neglect the $-term when: 



dH 



\V' -^V"\ < (3 + 2u) (dH'^^ , 
^ {2V - <i>V') < (3 + 200) (9H^) . 



(31) 



These conditions may be used to solve for when the slow-rollover approximation breaks down; 
that is, solved for values of $ for which the left-hand side of each inequality roughly equals the 
right-hand side (rather than being much less than it). Since $ = ^^/(8a;), the potential for our 
particular model may be written F($) = A/4 (8a;$ — v'^)^, which leads to: 



\V' -^V"\ ^AXujv'^. 



(32) 



After re-expressing H in terms of $ (see eq. (5)), we find the value of the field (<1>m) for which the 
consistency of the slow-roll approximation breaks down to be: 



V 

8u 



1 + 



37 



1 + 



1- yrre^ 
37 



1/2 



(33) 



If a; = 500, (^^d = 0.98?;. It is interesting to compare this with the result for based on the second 
condition for slow-rollover. This condition leads to the assignment: 

/72 1" 



%d = V 



1 + 



288 7 



1/2 



(34) 



In other words, the second condition doesn't break down until cf) > v\ Thus, the result based on 
the first condition will be used. 
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We may calculate a maximum reheating temperature for the model by finding the value of 
V{(j)) at the point where the field begins its damped oscillations around the true minimum of the 
potential. Taking this point to be (pbd leads to: 

viM = + 1 - yr+e^) • (35) 

If w = 500, V{(j)bd) = (3.2 X 10-^)Xv^, so TRH,max = V^^^iM = (0.13)Ai/^i;. Furthermore, if 
A ~ 10-^ then Trh max — (4.1 X 10 Assuming the simplest case, that cj) does not evolve after 
the end of inflation, then v = yJuj/2iTMp, which (for u) = 500) leads to T^ih, max ~ 4.5 x 10^'' GeV. 
We will consider possible interpretations of inflation at this energy scale below. 



4 Conclusions 



Induced-gravity Inflation, which combines properties from the "new inflation" schemes of 1982 |17| 
(such as a slowly rolling field leading to a second order phase transition) with characteristics from 
the original version of extended inflation (including a non-minimal R coupling), can lead 
to successful inflation with potentially acceptable limits on "fine-tuning". The lagrangian of eq. 
(2) requires only a qualitatively simple scalar potential associated with a single curvature-coupled 
scalar field; there is no need for adding special phenomenologically-inspired 'extra' terms by hand 
to £, as in [^, Induced-gravity Inflation can also get all of the 'work' of inflation done with only 
one boson, thereby helping to slow the proliferation of "specialty" bosons, each of which is invented 
to complete specific and unrelated tasks in the early universe. In addition to this simplicity, the 
model retains many of the motivations for extended inflation, based on the appearance of Brans- 
Dicke-like couplings in the low-energy effective theories for various Kaluza-Klein, superstring, and 
supergravity theories (see ||l8[|). 
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The calculation of A in this paper depends only on the amplitude of density perturbations, Sh- 
Yet the character of the spectrum of perturbations may also help to rule out various inflationary 



schemes 13, 21]. Ordinary extended inflation, for example, predicts a rather steep tilt away 
from a scale-invariant (Harrison-Zel'dovich) spectrum of density perturbations, which appears to 
contradict COBE data. Determination of the tensor mode contributions versus scalar modes in 
the density perturbation spectrum of Induced-gravity Inflation is made more difficult because of 
the deviation of eq. (7) from a simple power-law solution, and is the subject of further study. At 
early times at least, when the evolution of a(t) is roughly proportional to t"^, the present model 
would yield a tiny tilt away from scale- invariance: when a{t) oc V, the spectral index goes as 
n = 1 — 2/(7 — 1) [|lO|, which in this case (with lo ~ 500) would give n = 0.996. For more on 
the possibility of "observing" the inflaton potential based on the contributions from tensor mode 
perturbations, see |l^, 

One point of concern for Induced-gravity Inflation is the scale at which it operates: unlike most 
other inflationary schemes, which study phase transitions associated with the breaking of a GUT 
symmetry (at an energy of around 10^^ to 10^^ GeV), Induced-gravity Inflation is associated with 
the Planck scale. This could lead to conffict with the value of V{(p) at the time of last horizon- 
crossing. Refs. [0, m show that present COBE data appear to limit V^/^{U) ~ (3 - 4) x 10^^ 



GeV, which, for the present model (with lv ~ 500 and A ~ 10 ), would require 0(t*) to be very 



close to V. Yet, as pointed out in [13|, uncertainties in the data lead to an entire order of magnitude 



range in the value for y^/^(t*), so the present model cannot be ruled out by these COBE results. 



For more on constraints on the energy scale of inflation, see |19|. 

A theoretical difficulty for Induced-gravity Inflation stemming from its high energy scale is how 
to combine it with a "realistic" particle physics sector. (Recent work with extended technicolor as 
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a means of achieving Planck-scale unification of gauge couplings |22| might offer a means of con- 
necting an Induced-gravity Inflation model with realizable particle physics models.) Yet what it 
might lose on the particle side, it gains on the gravitational side: it should be much easier to relate 
the present model to a specific higher-energy gravitational theory. Or the model might be useful 
as part of a "Double-Inflation" scheme, in which the Induced-gravity phase transition (which, as 
we have seen above, could solve the flatness and horizon problems rather easily, and lead to an 
acceptable amplitude of density perturbations) is followed by a related GUT transition at a lower 
energy (which would then only need to solve the monopole density problem, so the requirements for 
this second epoch of inflation would be greatly relaxed). (For earlier attempts to use the original 
model of extended inflation in a double- inflation scenario, see pO[| .) Although these details have yet 
to be worked out, the prospect of a well-motivated inflationary scenario which requires A ~ 10^^ 
rather than ~ 10^^^ remains an encouraging result. 
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Figure Captions 



1. Plot of A as a function oi fi = (rrip/ME), based on eq. (26) with cu = 500. The vertical 
scale is in units of 10^^, and the assumption that aend/a(i*) ~ has been used. The intercept 
at At = is A = 3.8 X 10"^. 

2. Plot of A(a;) versus uj, based on eq. (26) with /x — 0. The vertical scale is in units of 10~^, and 
the assumption that aend/0'{t*) — has been used. Note that A reaches its maximum value of 
5.4 X 10-6 at cjcr = 240; the value at a; = 500 is A(500) = 3.8 x lO^^. 

3. Plot of A versus a, where acnd/a{t*) = e°, based on eq. (26) with /x — > and oj = 500. The 
vertical scale is in units of 10"^. 
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